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System and Methodology for Adaptive, Linear Model 
Predictive Control Based on Rigorous, Nonlinear Process 

Model 

1 . Field of the Invention 

This invention relates to the control of a 
multivariable, time-varying, industrial process that 
exhibits either linear or nonlinear response 
characteristics and more particularly the on-line 
execution of an adaptive, linear model predictive 
10 controller derived from a rigorous, nonlinear model of 
such processes. 

2 . Description of the Prior Art 

There are many multivariable, time-varying, 
industrial processes that exhibit either linear or 
nonlinear response characteristics. Several examples of 
such processes are a distillation column, a separator 
train, a catalytic cracking unit, a chemical reactor, and 
a utility boiler. Control of these processes using 
traditional multivariable control techniques with fixed 

20 model representations is difficult because of the variable 
process conditions that result from either process 
nonlinearities, measured or unmeasured disturbances, 
feedstock changes, or operator-induced changes. 

The standard practice for advanced industrial 
process control in processes of the type described above 
is to use linear, multivariable, model predictive 
controller (MPC) software. See for example, the 

Setpoint, Inc. product literature dated- 1993 entitled 
"SMC-Idcom: A State-of-the-Art Multivariable Predictive 

30 Controller"; the DMC Corp. product literature dated 1994 
entitled "DMC™: Technology Overview"; the Honeywell Inc. 
product literature dated 1995 entitled "RMPCT Concepts 
Reference"; and Garcia, C.E. and Morshedi, A.M. (1986), 
"Quadratic Programming Solution of Dynamic Matrix Control 
(QDMC)", Chem. Eng. Commun. 46: 73-87. The typical MPC 
software allows for model scheduling (i.e. changing the 
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model gains and/or dynamics) to improve control 
performance when operating on a nonlinear and/or time- 
varying process. The controller uses new models that are 
generally calculated in an off-line mode, or may be 
calculated by an adaptive algorithm that uses recent 
operating data. 

If the models are calculated off-line, then the 
controller requires additional on-line logic to determine 
which set of model parameters should be used at the 

10 current time. This logic is often difficult to develop 
since it may depend on numerous operating variables. For 
large problems there may be a significant number of 
different models required in order to improve 
performance. Furthermore, there is a possibility that 
invoking a particular set of model parameters for certain 
operating points will eventually lead to unstable 
operation. Thus, the off-line model identification task 
is extremely time-consuming and expensive and its 
implementation on-line is not proven robust over wide 

20 operating ranges. Therefore, in most applications a 
simplified set of models is defined (e.g. high feed rate 
and low feed rate models; or winter, and summer operation 
models) . 

The control performance using a simplified set of 
models is not much better than what can be achieved using 
a single model. In fact, controllers developed using 
this methodology do not usually remain on-line for 
processes demonstrating a high degree of non-linearity. 

Adaptive model identification is another way to 
30 modify linear models to capture current operating 
conditions. This procedure uses recent operating data to 
automatically adjust the models on-line. One of the 
biggest problems of this approach is inherent from the 
dual control principle, which essentially states that the 
uncertainty in model identification increases as the 
control performance improves (i.e. as the control 
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uncertainty decreases) . Thus, it is always much more 
difficult to obtain accurate new models when the 
controller is running since the input signals are 
correlated and the output signals have small deviations. 
In fact, the identification process often fails because 
of the loss of persistent excitation in the input 
signals. See Astrom, K.J. and Wittenmark, B. (1995) 
"Adaptive Control", Addison-Wesley Publishing Company, 
Inc., pp. 69-70 and 473-477, for a detailed discussion of 

10 this problem. This inherent mathematical problem limits 
the success of adaptive, linear model identification. 
Current algorithms use various methods to turn off or 
heavily filter the parameter changes in the adaptive 
identification. Otherwise, poor models could be selected 
because of the high model uncertainty and ultimately the 
control performance would suffer. 

There are some types of industrial processes that 
cause additional problems for empirical model 
identification methods. For example, a batch process 

20 could operate with multiple product runs during which the 
model gains change by orders of magnitude (e.g. 
polyethylene, polypropylene processing) , yet the batches 
may not reach equilibrium. Linear, empirical methods 
cannot obtain accurate models for such processes since 
the necessary operational data is not available. 

Model identification problems also occur for 
processes that do not allow plant testing near dangerous 
or unstable operating conditions (e.g. catalytic 
cracking, reactors) . Controllers for such processes try 

30 to prevent violation of constraint limits that could lead 
to unstable behavior, however the benefits often increase 
by operating closer to the limits. Therefore to improve 
performance and maximize benefits, the controller needs 
to have multiple models available for all the operating 
regions. This is not possible with linear, empirical 
methods. A rigorous, nonlinear process simulation model 
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could, however, provide these multiple models for the 
controller . 

The same type of model identification problems can 
occur for non-dangerous processes that have product 
specification limits (e.g. distillation columns). The 
plant manager does not want to generate off-spec product, 
so the testing must be conducted such that the purity 
remains within acceptable limits. However, greater 
benefits are typically obtained by operating closer to 

10 these limits. Thus, the empirical methods cannot 

generate accurate models for operating points that are 
close to or beyond the limits. High-purity separation 
processes present the most significant modeling problems 
as very nonlinear behavior may occur over the different 
operating regions . 

Another disadvantage of empirical MPC identification 
methods occurs when the controller is used together with 
real-time optimization (RTO) since the models are not 
consistent. RTO generally uses a more rigorous, nonlinear 

20 model while MPC is either using a fixed linear model 
obtained from plant testing, or an adaptive, linear model 
obtained from on-line data. In either case, the MPC 
models are not derived from the RTO model (nor vice 
versa) , so this inconsistency often leads to significant 
performance degradation. See for example the article by 
Y.Z. Friedman entitled "What's wrong with unit closed 
loop optimization?" which appeared in the October 1995 
issue of Hydrocarbon Processing- at pp. 107-116. 

Another means of treating nonlinear, time-varying 

30 processes is to use robust control methods. The idea 
behind robust control is to use a single process model, 
but to tune the controller by accounting for modeling 
uncertainties. The control design is tested against a 
range of expected operating conditions and is retuned 
until the performance is fairly consistent (or robust) 
over all conditions. The robust tuning method reduces the 
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sensitivity of the controller to model error. Honeywell 
has applied this technique to its RMPCT (Robust 
Multivariable Predictive Control Technology) product as 
described in the product literature referenced above. 

However, there are some drawbacks to using robust 
control techniques. In particular, by restricting the 
controller to a single linear model, robust tuning often 
results in slow, sluggish control performance over the 
typical operating range because it attempts to improve 
the worst-case performance. Thus, the robust tuning often 
degrades performance during periods of small modeling 
error. This occurs because the expected range of modeling 
uncertainties is often too large for a single model and 
single set of tuning parameters. It is possible to vary 
the tuning, but this becomes a difficult on-line 
implementation task. The real solution is to allow 
multiple models for processes with significant 
nonlinearities and/or time-varying behavior. 

The patent prior art includes many examples of the 
use of model-based control systems employing both linear 
and nonlinear methodologies for control. Most of the 
prior art in MPC refers to linear controllers [see for 
example, U.S Patent Nos. 4,349,869 to Prett et al.; 
4,616,308 to Morshedi et al.; 5,351,184 to Lu et al.; and 
5,572,420 to Lu] . To handle nonlinear, time-varying 
processes, these controllers may use gain scheduling, 
adaptive model estimation, or robust control tuning. 
These approaches typically encounter implementation 
problems and/or performance degradation for the types of 
processes and operating conditions described previously. 

There have been a few patents issued for nonlinear 
model-based control methodologies. In particular, U.S. 
Pat. No. 5,260,865 to Beauford et al. describes a 
nonlinear model-based control strategy for a distillation 
process which employs a nonlinear model to compute liquid 
and vapor flow rates required for composition control. 
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Sanner and Slotine (US Pat. No. 5,268,834) employ a 
neural network together with other nonlinear control 
strategies to provide adaptive control of a plant. 
Bartusiak and Fontaine (U.S. Pat. No. 5,682,309) developed 
a reference synthesis technique in which the controller 
attempts to make a nonlinear plant follow a specified 
reference trajectory. US Pat. No. 5,740,033 to Wassick et 
al. describes an MPC that employs a real-time executive 
sequencer and an interactive modeler to find the optimized 

10 set of control changes for a nonlinear process. 

Large, nonlinear control problems are difficult to 
solve in an on-line operating environment. The solver 
must be fast and robust. It is unknown whether these 
inventions have been able to achieve successful 
implementations of on-line, nonlinear control. 

Summary of the Invention 
The invention is a methodology for process modeling 
and control and the software system implementation of this 
methodology, which includes a rigorous, nonlinear process 

20 simulation model, the generation of appropriate linear 
models derived from the rigorous model, and an adaptive, 
linear model predictive controller that utilizes the 
derived linear models. A state space, multivariable, model 
predictive controller (MPC) is the preferred choice for 
the MPC since the nonlinear simulation model is 
analytically translated into a set of linear state 
equations. Thus, a state space MPC simplifies the 
translation of the linearized simulation equations to the 
modeling format required by the controller. However, the 

30 state space methodology also allows the calculation of 
various other MPC modeling forms such as transfer 
functions, impulse response coefficients, and step 
response coefficients. Therefore, the methodology is very 
general in that any model predictive controller using one 
of the above modeling forms can be used as the controller. 
The methodology also includes various modules that 
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improve reliability and performance. For example, there 
is a data pretreatment module used to pre-process the 
plant measurements for gross error detection. A data 
reconciliation and parameter estimation module is then 
used to correct for instrumentation errors and to adjust 
model parameters based on current operating conditions. 
The full-order state space model can be reduced by the 
order reduction module to obtain fewer states for the 
controller model. Automated MPC tuning is also provided to 
improve control performance. 

The present invention is embodied as a method for 
controlling a process. The method has the steps of: 

a) \ receiving plant measurement variables from a 
regulator^ control system; 

b) Applying the plant measurement variables to 
update one o\ more variables of a nonlinear model; 

c) linearizing the updated nonlinear model; and 

d) passing a MPC format model converted from the 
linearized model t\ a model predictive controller. 

The invention i^s also embodied as a method for 
controlling a process tfc^at has the steps of: 

a) converting a not^linear model to a linear model 
for operation of a model predictive controller; 

b) using the nonlinear model in a real time 
optimizer to compute targets Bor the model predictive 
controller; and 

g-) — —p ass i ng tE^ rerarget s Lo Lhe cofrtrol3 re-fr. 
T^he invention- is further embodied as a method for 
control Irfrg^a process that has the steps of: 

a) conv^s4^ing at least one submodel of a nonlinear 
model having two oKmore submodels to a linear model, 
each of the two or more^u^models having a predetermined 
one of two or more modeT^oredictive controllers 
associated therewith, the linear mod^l for operating the 
associated one of the two or more controllers; 

b) using the nonlinear model in aX.real time 
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optimizer to compute targets for all of the two or more 
model predictive controllers , a predetermined siaifget of 
the computed targets associated with a re^p€ctive one of 
the two or more controllers; ^^^^ 

c) passing each of tja^predetermined subsets of 
the computed targets ^^s^ociated with a respective one of 
the two or mope^model predictive controllers to the 
associated^oite of the two or more controllers; and 

&r passing the linear model to the associated one 
o^the two or more controllers. 

The invention is also further embodied as a method 
for controlling a process that has the steps of: 

a) Applying simulation stimuli to update one or 
rnOgre variables of a nonlinear model; 
linearizing the updated nonlinear model; 
converging the linearized model to a MPC format 
model; abnd 

passing blje MPC format model to a model 
predictive cWitroller. 

invention is \embodied as an apparatus for 
controlling a process \hyaving process measurement 
variables associated therewith, \he apparatus comprising: 
a) a digital procesy 

a model predictivp^cof^roller having a model 
for the process"^! 1 y r e-i 

a simulation environment \ routine having a 
nonlinear modeb. I therein, \ the simulation 
environment routine executed 
processor for: 

(i) applying process measuremeVb variables 
to update one or more variables of the 
nonlinear model; 

(ii) linearizing the updated nonlineai\model; 
and 

(iii) passing a MPC format model converted 
from the linearized nonlinear model 



b) 
c) 

d) 

The 



b) 
c) 



>y the digital 
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the model predictive controller, 
v The invention is also embodied as an apparatus for 
controlling a process having process measurement 
variables associated therewith, the apparatus comprising: 

a) axiigital processor; 

b) a model predictive controller having a model for 
the^process therein; and 

c) a simulation environment routine having a 
nonlineax model therein, the simulation 
environment /x?out\ine executed by the digital 
processor foi 



(i) 
(ii) 

(iii) 



con 
nonl 
usin 
opt 



ng J 
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the linear model to the 
_p_n|ear fltodel; and 

the norridnear model in a real time 
.zer to cbjnpute targets for the 
model predictive coj^troller; and 
passing the target^ to the model 
predictive controller . 
The invention is further embodied as an apparatus 
for controlling a process having process measurement 
variables associated therewith, the apparatus comprising: 

a) a digital processor; 

b) two or more model predictive controllers each 

having an associated submodel of a linear 
model having two or more submodels for the 
process therein; and 

c) a simulation environment routine having the 

nonlinear model therein, the simulation 
environment routine executed by the digital 
processor for: 

(i) converting at least one submodel of a 
nonlinear model having two or more 
submodels to the associated linear 
submodel; 

(ii) using the nonlinear model in a real time 
optimizer to compute targets for all of 
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the two or more model predictive 
controllers, a predetermined subset of 
the computed targets associated with a 
respective one of the two or more model 
predictive controllers; 

(iii) passing each of the predetermined 
subsets of the computed targets 
associated with a respective one of the 
two or more model predictive controllers 
to the associated one of the two or more 
controllers; and 

(iv) passing the linear submodel to the 
associated one of the two or more 
controllers . 

invention is also further embodied as an 
apparatiiks for controlling a process having process 
measurement variables associated therewith, the apparatus 
comprising: 

a) a digi*&al processor; 

b) a model predictive controller having a model for 

the process therein; and 

c) a simulation x^j environment routine having a 
nonlinear mocfei^ i\ therein, the simulation 
environment routiii^c executed by the digital 

processor for: 

(i) applying slimulata>on stimuli to update 
one or more variables of the nonlinear 
model; 

(ii) linearizing the updated \onlinear model; 
and 

(iii) passing a MPC format mod£\L converted 
from the linearized nonlineaK model to 
the model predict ivo -^-con'tToTll 

i 

Description of the Drawing 
Fig. 1 shows in a block: diagram form the various 
software components of the invention in the multiple 
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operating environments and the data flow between 
components . 

Fig. 2 shows a variation of the embodiment of Fig. 1 
in which the model change detection module is part of the 
Model and Tuning Environment. 

Fig. 3 shows a variation of the embodiment of Fig. 1 
using a different type of MPC controller, i.e. not a state 
space MPC controller. 

Fig. 4 shows an extension of the embodiment of Fig. 1 
that includes the described methodology together with a 
real-time optimization module in the Simulation and 
Optimization Environment. 

Fig. 5 shows an extension of the embodiment of Fig. 1 
similar to Fig. 4 with a real-time optimization module 
that is setting targets for multiple MPCs, so there are 
multiple MPC Environments and multiple Model and Tuning 
Environments . 

Fig. 6 shows a variation of the embodiment of Fig. 1 
used for off-line design of an adaptive, linear MPC based 
on a rigorous, nonlinear model. 

Description of the Preferred Embodiment (s) 

The embodiment of the complete computer software 
system 50 shown in Fig. 1 includes a number of different 
components operating in multiple environments. These 
multiple environments include a Simulation Environment 10, 
which receives inputs from a Plant Environment 20; and a 
Model and Tuning Environment 30 which receives inputs from 
the Simulation and Plant Environments 10, 20 and from a 
MPC Environment 40. The MPC Environment 40 receives 
inputs from the Model and Tuning Environment 30 and from 
the Plant Environment 20 and provides outputs to the Plant 
Environment 20 and to the Model and Tuning Environment 30. 

The Plant Environment 20 includes regulatory control 
22 to control the operation of process 24. As is well 
known regulatory control 22 typically comprises a process 
automation system that includes associated instrumentation 
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such as sensors for measuring temperature, pressure and 
flow; and analyzers. 

Each of the environments 10, 30 and 40 include the 
following multiple software components: 
Simulation Environment 10 

• Nonlinear simulation model 12 

• Data Pretreatment 14 

• Data reconciliation and parameter estimation 16 

• Model change detection 18 

10 • Linear model export 19 

Model and Tuning Environment 30 

• State Space Conversion 32 

• Order reduction 34 

• Discrete-time conversion 35 

• Model Evaluation 36 

• Automated MPC tuning 38 
MPC Environment 40 

• Model transfer smoothing 4 2 

• State estimator and controller 44 

20 The Simulation Environment 10 includes five 

components or modules: nonlinear model 12, data 
pretreatment 14, data reconciliation and parameter 
estimation 16, model change detection 18, and linear model 
export 19. The nonlinear model 12 is a representation of 
the industrial process as a set of mathematical equations 
(ordinary differential, partial differential, integral, 
and/or algebraic equations) . The nonlinear model 

typically has the following general form: 

g (x(t)Mt),y(0,0) = o 

30 where 

/ represents time, 

x(t) is the state vector, 
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x(t) is the time derivative of the state vector, 
u(t) is the input vector, 
y(t) is the algebraic vector, 
0 is the parameter vector, 
/ defines the differential equations, and 
g defines the algebraic equations. 

The plant data from the Plant Environment 20 is 
processed by the data pretreatment module 14, which is 
designed to detect and remove gross errors from the plant 

10 data. The nonlinear model 12 and the data reconciliation 
and parameter estimation module 16 each receive the 
processed plant data from module 14. 

The data reconciliation and parameter estimation 
module 16 balances correcting sensor-based measurement 
errors against readjusting unmeasured disturbances (e.g. 
equipment attributes or feed characteristics) to provide a 
real-time estimate of the predicted plant behavior. This 
module receives measurements preprocessed (i.e. gross 
error removal) by module 14 from the plant instrumentation 

20 and control system (regulatory control 22 in the Plant 
Environment 20) and processes this information to create a 
better representation of the plant's performance. The data 
reconciliation procedure is required to correct for 
instrumentation errors, such as thermocouple drift, 
improper flowmeter calibration, etc. by altering the 
preprocessed measurement data. The parameter estimation 
procedure uses the reconciled measurement data to readjust 
model parameters (9) resulting from changes in unmeasured 
disturbances that impact model performance. The nonlinear 

30 model 12 then uses the reconciled data and new parameter 
estimates to estimate the values of the other plant 
variables . 

The model parameters may represent physical 
properties or special properties defined for modeling 
purposes (e.g. heat exchange coefficients, tank mixing 
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factors, aggregate feed composition properties, etc.). 
Since the parameters are closely associated with the 
actual process and plant equipment (as opposed to transfer 
function coefficients and step response coefficients) and 
since the data is also reconciled by module 16, this type 
of parameter estimation procedure generally provides more 
accurate on-line estimates than an adaptive, linear model 
identification method. The new model parameters are then 
passed to the nonlinear model module 12. 

The model change detection module 18 evaluates the 
changes in the model variables, which are defined herein 

as parameters (0) and dynamic variables (x, y, u) , and 
determines when a new linear model should be generated. 
The model change detection logic can be based on various 
criteria (e.g. when a single variable changes by more than 
a predetermined threshold or if the sum of the absolute 
values of all variable changes exceeds a different 
threshold, etc.). When the change detection module 18 
determines, using one of the criteria described above, 
that there is a change in the model variables, it sets a 
flag and sends it to the nonlinear model module 12. The 
nonlinear module 12 then activates the linear model export 
module 19 to analytically generate a set of linearized 
model equations from the nonlinear model of module 12. 

The software of the linear model export module 19 
provides symbolic differentiation of the nonlinear model 
equations to create linearized equations in the following 
form: 
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This set of equations is then passed to the Model and 
Tuning Environment 30 in a sparse matrix format. 

The Model and Tuning Environment 30 includes five 
modules: state space conversion 32, order reduction 34, 
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discrete-time conversion 35, model evaluation 36, and 
automated MPC tuning 38. The state space conversion module 
32 receives the implicit linear model equations from the 
rigorous Simulation Environment 10. These equations are 
then translated by the state space conversion module 32 
into a state space matrix model based on the control 
structure specified in the current controller model, i.e. 
the set of controller inputs (manipulated variables and 
disturbance variables) and outputs (controlled variables) 
10 received by Model and Tuning Environment 30 from MPC 
Environment 40. The resulting state space model has the 
following standard, continuous-time form: 

x(t) = Ax(t) + Bu(t) 
y{t) = Cx{t) + Du(t) 

The order reduction module 34 is then used to obtain 
an approximately equivalent state space model, in terms of 
state observability and controllability, with a reduced 
number of states. This procedure is extremely important 
for some practical applications as it may reduce the 
number of states by one or more orders of magnitude. The 
20 controller design is also simplified by reducing the 
number of states. 

The standard technique for order reduction is to 
convert the state space model to a balanced realization 
and then truncate some of the states (commonly referred to 
as balance and truncate method) . This method involves 
finding a balanced realization in which the 
controllability and observability grammians are equal, and 
the states are sorted from the most controllable and 
observable, to the least. The states that are much less 
30 observable and much less controllable than the others can 
then be eliminated (based on a specified threshold) to 
form the reduced order model. Alternatively, other well 
known order reduction methods could be used in place of 
the standard balance and truncate method. 

It should be appreciated that the order reduction 



E005640 16 

module 34 can be bypassed if desired as is shown in Fig. 1 
by dotted line 33. For example, such bypassing may occur 
if the simulation model has small number of states; or to 
retain same simulation state variables in state space 
model . 

The reduced order model is then converted to a 
discrete-time state space model by module 35 based on the 
desired sampling period, T, for the MPC. The MPC format of 
the model is: 

x(k + l) = A d x{k) + B oul u(k) + B ff d{k) 
y(k) = C d x(k) + D 0Ul u(k) + D jr d(k) 

where 

k represents the time step, 

x(k) is the state vector, 

u(k) is the manipulated variable vector, 

d(k) is the disturbance variable vector, and 

y(k) is the controlled process variable. 

Alternatively, the order reduction methods may be 
applied to a full-order, discrete-time state space model. 
In this alternative procedure, which is not shown in Fig. 
1, the full-order, continuous-time state space model would 
first be converted to a discrete-time model, and then 
order reduction techniques would be applied to the 
discrete-time model. 

Following the creation of a reduced order, discrete- 
time state space model, the model evaluation procedure of 
module 36 is used to evaluate the time and frequency 
domain characteristics of the model. Various open-loop and 
closed-loop performance criteria are used to compare the 
new reduced order model to the new full-order model and to 
the existing controller model, e.g. 

(01) Compare predictions for recent operating data 
(starting at the time when the parameter or variable 
change was detected) . 

(02) Compare open-loop step responses. 
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(03) Compare open-loop frequency responses . 

(CI) Compare closed-loop step responses. 

(C2) Compare closed-loop frequency responses. 

The model evaluation logic can assume one of various 
forms. For example, compute a weighted sum of the open- 
loop performance criteria (01, 02, and 03) and compare to 
a specified open-loop limit. If the weighted open-loop 
performance measure is less than the open-loop limit, then 
the reduced order model is rejected and a new reduced 
10 order model is calculated by the order reduction module 
34. The open-loop performance measure of this new reduced 
order model is then evaluated. 

If the weighted open-loop performance measure exceeds 
the open-loop limit, then proceed to compute a weighted 
sum of the closed-loop performance criteria (CI, C2) using 
the existing MPC tuning specifications and compare this to 
a closed-loop limit. If the weighted closed-loop 
performance measure exceeds the closed-loop limit, then 
the new reduced order model is accepted and passed to the 
20 MPC Environment 40 (with the existing MPC tuning specs) . 

However, if the weighted closed-loop performance 
measure is less than the closed-loop limit, then the 
automated MPC tuning procedure in module 38 is invoked to 
define new controller tuning specifications for this 
model . 

The automated tuning module 38 tries to find new 
tuning specifications for the MPC state estimator and 
controller 44. Various tuning methods can be used. For 
example,, one method is to find new tuning specifications 
30 that improve the closed-loop performance measure of the 
new model while matching tuning speed factors of the 
current model. The tuning speed factors are calculated for 
each controlled variable based on the closed-loop speed of 
response relative to the open-loop speed of response using 
predetermined typical manipulated variable moves. As a 
first step, this module reads information from the MPC 
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Environment 40 regarding the current model and tuning 
specif ications, including the tuning speed factors. The 
automated tuning module 38 then performs an iterative 
search to find the tuning specifications for the new model 
that maximize the weighted closed-loop performance measure 
while maintaining the same tuning speed factors as in the 
existing model. 

Thus, in the Model and Tuning Environment 30, it is 
possible to have iterations between the order reduction 

10 34, model evaluation 36, and automated MPC tuning 38 
modules. When the iterations are complete, a new model and 
possibly new tuning specifications are passed to the MPC 
Environment 40. If the number of iterations reaches a 
predetermined limit, then the iterations cease without 
passing a new model to the MPC Environment 40. 

It is also possible to use a simplified procedure for 
on-line model evaluation in order to reduce the 
computation time (e.g. just use performance measure 01, or 
01 and CI) . In that case, the off-line model evaluation 

20 procedure may be more complicated (such as the one 
described above) and the results may be used to help 
define the simplified on-line procedure (e.g. find 
appropriate thresholds for 01 and CI) . 

Whenever a new model and tuning specifications are 
accepted by the model evaluation module 36, the MPC 
Environment 40 receives the new model and tuning 
specifications. This information can either be 

•automatically downloaded to the MPC 44, or presented to an 
operator /engineer for approval prior to download. After a 

30 new model is downloaded to the MPC 44, the model transfer 
smoothing module 42 is invoked to smooth the transfer 
between the two models 

In one embodiment of the present invention, the MPC 
controller was embodied by an AdvaControl 3dMPC available 
from ABB. This controller is described in Sweden patent 
application serial number 9903592-5, the disclosure of 
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which is expressly incorporated herein by reference. In 
such an embodiment, the model transfer smoothing module 
uses the observer or state estimator tuning specifications 
passed from the Model and Tuning Environment 30 to 
calculate a new gain vector for the state estimator of the 
new model and then uses past operating data to smooth its 
initialization. The MPC state estimator updates the 
estimates of the states based on linear model predictions 
and current operating data. 

10 In this embodiment of the present invention the MPC 

controller uses the state estimates and finds the optimal 
solution to a control problem of minimizing a weighted sum 
of manipulated variable moves plus controlled variable 
deviations from setpoints, and violations of constraints. 
The control problem also includes hard absolute and rate 
of change constraints for the manipulated variables. A 
move horizon is defined for each manipulated variable and 
a prediction horizon is defined for each controlled 
variable. The MPC solution defines the proposed moves for 

20 each manipulated variable during its move horizon. For 
each manipulated variable, the first calculated move is 
added to its current value and then sent to the plant 
environment to be implemented in the regulatory control 
module 22. The manipulated variable values represent new 
targets for the regulatory controllers (usually PID 
controller setpoints) or outputs to actuators. These 
regulatory control changes affect the process, then new 
data is- sent to each of the other three • environments 
(Simulation 10, Model and Tuning 30, and MPC 40), and the 

30 entire procedure is repeated. 

Figure 2 shows a variation 100 of the preferred 
embodiment of Figure 1 in which the model change 
detection module 18 of the Simulation Environment 10 of 
Fig. 1 is now included as module 31 in the Model and 
Tuning Environment 30' of Fig. 2. In this representation, 
the model change detection module 31 either receives 
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model predictions from the MPC Environment 40 and 
compares these to the measurement data from the Plant 
Environment 20 to determine the model prediction errors 
or depending on the type of MPC receives the prediction 
errors directly from the MPC. 

The model change detection module 31 then evaluates 
the model prediction errors and determines when a new 
linear model should be generated. The change detection 
logic of module 31 can be based on various criteria (e.g. 

10 when a single prediction error exceeds a predetermined 
threshold or the sum of the absolute values of all 
prediction errors exceeds a different threshold) . When 
the change detection logic determines that there is 
significant model error, it sets a flag and sends it to 
the Simulation Environment 10' . 

The linear model export module 19 of the Simulation 
Environment 10' then provides linearized equations of the 
nonlinear model of module 12 to the state space 
conversion module 32 of the Model and Tuning Environment 

20 30' . Other than this distinction, the remaining modules 
function as described in the preferred embodiment and 
have reference numerals which are identical to those used 
in Fig. 1. 

Figure 3 shows a variation 110 of the preferred 
embodiment of Figure 1 using a non-state space MPC 
wherein the environments and modules that are the same as 
those in Fig. 1 have reference numerals which are 
identical to those used in Fig. 1 for those environments 
and modules. The MPC Environment 40' of Fig. 3 includes' a 
30 model predictor and controller 44' instead of the state 
estimator and controller 4 4 shown in the embodiment of 
Fig. 1. The Model and Tuning Environment 30' has a 
conversion and model calculation module 34' that replaces 
the order reduction and discrete-time conversion modules 
34, 35 of the embodiment of Fig. 1. The conversion and 
model calculation module 34' includes a transformation of 
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the full order continuous-time state space model provided 
by the state space conversion module 32 to a discrete- 
time model, and the calculation of the model parameters 
in their corresponding MPC format. 

In a further variation, not shown in Fig. 3, the 
conversion and model calculation module 34 ' could 
transform a reduced order state space model provided by 
the order reduction module 34 of Fig. 1 to a discrete- 
time model (instead of converting the full order state 

10 space model) . Another variation, also not shown in Fig. 
3, is to add a module that reduces the complexity of the 
model generated by the conversion and model calculation 
module 34'. For example, a module could be added to 
reduce the order of a transfer function matrix 
representation of the model. 

A further variation, also not shown in Fig. 3, is to 
first reduce the complexity of the linear model and then 
convert it to the MPC format required to update the MPC 
model. For example, it is possible to only extract the 

20 linearized steady-state (DC) gains from the rigorous 
nonlinear simulation model and then apply only these 
gains (without changes in the dynamics) to the MPC model. 
If the new model results from any of the procedures using 
the various embodiments for the conversion and calculation 
module 34' described above in connection with Fig. 3, that 
model is referred to as a simplified MPC format model. 

Figure 4 shows an extension 120 of the preferred 
embodiment of Figure 1- to also include a real-time 
optimization module 17 in the modified Simulation and 

30 Optimization environment 10 f . The environments and 
modules in Fig. 4 that are the same as those in Fig. 1 
have reference numerals which are identical to those used 
in Fig. 1 for the corresponding environments and modules. 
The real-time optimization module 17 of Fig. 4 uses the 
rigorous, nonlinear process model data from module 12 
(the data could be read from a shared memory location) . 
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This model data is reconciled and includes updated model 
parameters from module 16. 

The real-time optimization module 17 then solves for 
the values of the optimization variables that optimize 
its objective function (e.g. to minimize operating 
costs) . The output data of the real-time optimization 
module 17 typically represents targets for the MPC 
controller (e.g. steady-state targets for some or all of 
the manipulated variables and setpoint targets for some 

10 or all of the controlled variables) . This data is passed 
to the state estimator and controller 44 of the MPC 
Environment 40 for implementation. The output data of the 
real-time optimization module 17 is also passed to the 
nonlinear model 12 within the Simulation and Optimization 
Environment 10'. Alternatively, the output of module 17 
could be passed to a memory location that is shared with 
the nonlinear model 12. 

Depending on the required and available hardware 
resources, the optimization extension of the preferred 

20 embodiment may include a separate Optimization 
Environment, not shown in Fig. 4, for the real-time 
optimization module 17. In this case, the Optimization 
Environment would communicate with the Simulation 
Environment to utilize the same rigorous, nonlinear model 
data. The real-time optimization module 17 would send its 
output to both the MPC Environment and to the Simulation 
Environment . 

Figure 5 shows a modified version 130 of Figure 4 in 
which the nonlinear model 12 and the real-time 
30 optimization module 17 of the Simulation and Optimization 
Environment 10 1 communicate with multiple MPC 
Environments 40 1 and multiple Model and Tuning 
Environments 30 ' . In this embodiment, there is a single 
nonlinear model 12 representing the complete plant and a 
single real-time optimization module 17. However, there 
can be more than one MPC, each MPC associated with a 
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section of the complete plant model (i.e. a submodel). 
Each MPC uses a model representing the corresponding 
plant submodel to control this portion of the plant. 
Thus, if there are N MPCs, then as is shown in Fig. 5, 
the MPC Environment 40 1 comprises 1 to N such 
environments and the Model and Tuning Environment 30 1 
also comprises 1 to N such environments, each of the 1 to 
N Environments 30' , 40' associated with a respective one 
of the N MPCs. The optimization module 17 provides 
10 targets for some or all of the N MPC controllers. 

The model change detection module 18' of the 
Simulation and Optimization Environment 10 1 determines 
13 when there is a significant change in one or more of the 

|y N submodels. When this occurs, the model change detection 

: p module 18 1 sets flags for each of the N submodels in 

which there are significant changes. The linear model 

: i [ 

jjg export module 19 1 then provides new linearized submodel 

P eguations to each of the N Model and Tuning Environments 

]^ 30' whose corresponding change detection flag is set. 

:p 20 In another variation not shown in Figure 5, there 

~ could, as is shown in Fig. 2, be a model change detection 

q module 31 1 in each Model and Tuning Environment 30 f 

Q instead of having a single model change detection module 

in the Simulation and Optimization Environment 10 In 
this variation, each model change detection module 
determines when to generate a new linear submodel for its 
corresponding MPC. 

■Figure 6 shows a variation 160 of the embodiment of 
Fig. 1 used for off-line design of an adaptive, linear MPC 
30 based on a rigorous, nonlinear model. In this embodiment, 
the Plant Environment 20 of Fig. 1 is replaced by a 
Stimulus Environment 60 which provides the simulated 
operator actions and process disturbances. The term 
M simulation stimuli" is used herein to refer to the 
simulated operator actions and process disturbances. The 
simulation stimuli can be generated automatically from a 
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pre-determined event file or randomly from an event 
generator, or spontaneously from the user or can be 
generated from any combination of the foregoing sources. 
The simulation stimuli is sent to the nonlinear model 12 
of the Simulation Environment 10. The nonlinear model 12 
provides simulated plant measurement data for the data 
pretreatment module 14 and for the Model and Tuning 
Environment 30 and the MPC Environment 40. The MPC 
Controller 4 4 provides targets for the regulatory 
controllers which are simulated within the nonlinear model 
12. The remaining modules and data transfers between 
modules function in exactly the same manner as in the 
embodiment of Fig . 1 . 

It is to be understood that the description of the 
preferred embodiment ( s ) is (are) intended to be only 
illustrative, rather than exhaustive, of the present 
invention. Those of ordinary skill will be able to make 
certain additions, deletions, and/or modifications to the 
embodiment (s) of the disclosed subject matter without 
departing from the spirit of the invention or its scope, 
as defined by the appended claims. 



